Accepted for publication in ApJL 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



GALACTIC EMISSION AT 19 GHZ 

Angelica de Oliveira-Costa 1,2,3 , Max Tegmark 3 ' 4 , Lyman A. Page 1 & Stephen P. Boughn 5 

Accepted for publication in ApJL 

ABSTRACT 

We cross-correlate a 19 GHz full sky Cosmic Microwave Background (CMB) survey with other maps 
to quantify the foreground contribution. Correlations are detected with the Diffuse Infrared Background 
Experiment (DIRBE) 240, 140 and 100/iin maps at high latitudes (|6|>30°), and marginal correlations 
are detected with the Haslam 408 MHz and the Reich & Reich 1420 MHz synchrotron maps. The 
former agree well with extrapolations from higher frequencies probed by the COBE DMR and Saskatoon 
experiments and are consistent with both free-free and rotating dust grain emission. 

Subject headings: cosmic microwave background - methods: data analysis 



1. INTRODUCTION 

One of the major challenges in any Cosmic Microwave 
Background (CMB) anisotropy analysis is to determine 
, the fraction of the observed signal due to diffuse Galac- 
tic emission. Three components of Galactic emission have 
been firmly identified: synchrotron and free-free radiation, 
which are important mainly at frequencies below 60 GHz, 
and thermal emission from dust particles, which is impor- 
tant mainly at frequencies above 60 GHz (see, e.g., Weiss 
■ 1980; Bennett et al. 1992; Brandt et al. 1994; Tegmark 
and Efstathiou 1996). In principle, these three compo- 
nents can be discriminated by their frequency dependence 
and morphology. In practice, however, there is no com- 
ponent for which both the frequency dependence and the 
spatial template are currently well known (see, e.g., Kogut 
et al. 1996a, hereafter K96a, and references therein). 

The cross-correlation technique allows one to estimate 
the rms level of Galactic emission present in a CMB 
map. For instance, an analysis of high latitude Galactic 
emission in the COBE DMR map gave rms estimates of 
(3.4±3.7) for synchrotron and (2.7±1.3) fiK for dust 
emission at 53 GHz on a 7° scale (Kogut et al. 1996b, 
hereafter K96b). A third component, correlated with the 
DIRBE maps but decreasing with frequency, was detected 
at the level of (7.1±1.7) /jK and tentatively identified as 
free- free emission (K96b). Although this component was 
also detected in the Saskatoon maps (de Oliveira-Costa 
et al. 1997, hereafter dOC97) on a 1° scale, at a level 
(17.5±9.5) [iK at 40 GHz, it is still not clear if it is due 
to free-free emission. An analysis of Owens Valley Ra- 
dio Observatory (OVRO) data at 14.5 GHz showed such a 
component at the level of 203 /iK in a small sky region on 
scales of 7'-22' (Leitch et al. 1997, hereafter L97), but this 
level is substantially higher than free-free emission esti- 
mates based on Ha images (Gaustad et al. 1996; Simonctti 
et al. 1996). Suggested explanations include the presence 
of high temperature gas (L97) and rotating dust grains 
(Draine & Lazarian 1998, hereafter DL98), but more data 
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Fig. 1. — 19 GHz survey (top), 408 MHz Haslam synchrotron 
template (middle), and DIRBE 100/^m dust template (bottom). 

is needed to settle this issue. The purpose of this letter is 
to evaluate the Galactic contribution to the 19 GHz full 
sky map (see Figure [|) by cross-correlating it with the 
DIRBE dust maps, and with the Haslam and Reich & Re- 
ich synchrotron maps. 

2. METHOD 

The 19 GHz map consists of N = 24576 pixels with 
sky temperatures j/j and noise Uj. We assume that this 
map is a superposition of CMB fluctuations and a Galac- 
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tic component whose angular distribution is traced in part 
by an external data set. Representing these contributions 
by jV-dimensional vectors, 



y = n + X C MB + aX G al + YGal , 



(1) 



where x GMB is the contribution of the fluctuating com- 
ponent of the CMB, x Gal is the brightness fluctuations of 
the Galactic template map (not necessarily in temperature 
units), a is the coefficient that converts units of the Galac- 
tic template into antenna temperature and y Gal represents 
any residual Galactic contribution which is uncorrelated 
with x Gal . We consider n and xcmb to be random vari- 
ables with zero mean, i.e. (xcmb) = (n) = 0, and X-Gai 
and yGai to be constant vectors. Thus the data covariance 
matrix is given by 

C ee (yy T ) - <y)(v T ) = (x CMB x T GMB ) + <nn T ), (2) 

where (xcmb^-cmb) * s the covariance matrix of the CMB 
and (nn T ) is the noise covariance matrix. The noise in 
the 19 GHz map is approximately uncorrelated and has 
an rms amplitude of (Tj~2 mK. Therefore, the covariance 
matrix of this map is 



C i:j w (mrij) 



rfS, 



Minimizing \ 2 = (y - ax Gai ) T C x (y 
the minimum-variance estimate of a, i.e. 



^Gal 



^Gal 



(3) 

ax Ga/ ) yields 



(4) 



with variance 



A« 2 



Gal 



XGal) 
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Note that unlike the case in dOC97, there is little contri- 
bution from chance alignments between the CMB and the 
various template maps, since the CMB contribution to C 
is negligible. If the the noise is correlated or the anisotropy 
signal is significant, then our C used above will differ from 
the true covariance matrix, denoted C Although eq. ^ 
still provides a reasonable and unbiased estimate of a 1 , its 
variance will be larger than implied by eq. (j^), given by 



Aa 2 



(xg n , C- 1 C C 1 XGal) 



( X Ga2 



C- 1 x Gai ) 2 



(6) 



In the next section, an estimate of C is made from the 
data and the corresponding variance of a is evaluated. 

3. DATA ANALYSIS AND RESULTS 

The 19 GHz map has an angular resolution of 3° FWHM 
and is stored in 1.3° x 1.3° pixels (Cottingham 1987; 
Boughn et al. 1992). The template maps are convolved 
with a 3° Gaussian beam and regions within 20° and 30° of 
the Galactic plane are excluded. To avoid contamination 
by zodiacal dust emission, data within 10° of the Eclip- 
tic plane are also excluded from the analysis; although, 
the results are found to be independent of this cut. Off 

1 Note that the diagonal approximation of eq. (0) reduces eq. (W) to 



the Galactic plane, the 19 GHz map is dominated by the 
CMB dipole (^=1). On the other hand, because of its pla- 
nar structure, emission associated with the Galaxy has 
a strong quadrupole {£=2) component. Therefore, the 
monopole, dipole, and quadrupole moments are removed 
from both the 19 GHz and template maps. As a conse- 
quence, the computed a depends only on correlated struc- 
ture in the maps on angular scales < 90°. 

3.1. Correlations and their Error Bars 

The 19 GHz map was cross-correlated with five differ- 
ent templates: two for synchrotron emission, the 408 MHz 
(Haslam et al. 1981) and 1420 MHz (Reich and Reich 1988) 
surveys; and three to study dust and free- free emission, the 
100, 140 and 240/im Diffuse Infrared Background Exper- 
iment (DIRBE) sky maps (Boggess et al. 1992). Table 1 
lists the coefficients d derived from eq. (|]) with errors com- 
puted from eq. (^). All three DIRBE templates show sig- 
nificant correlations with the 19 GHz map, while the two 
synchrotron templates are found to be only marginally cor- 
related. 

Also listed in Table 1 are the implied fluctuations in an- 
tenna temperature in the 19 GHz map, i.e. ST = aacai, 
where (7 Ga z is the rms of the template map. If a template 
map includes a distinct component that is uncorrelated 
with the 19 GHz map, then 5T is underestimated by a 
factor o-' Gal /aGah where o' Gal is the rms of the correlated 
component of the template map. For this reason, the <5T"s 
in Table 1 should be considered lower limits. 

As mentioned above, the errors listed in Table 1 were 
computed from eq. (^|) and are therefore lower limits to 
the error. We now describe a series of tests, performed to 
estimate the uncertainty in d due to correlated noise and 
other systematics. 

3.1.1. The Noise Correlation Function 

If the true noise correlation matrix is isotropic, C'^ — 
(niUj) — R(6ij)<JiO~j where 9ij is the angle between pixels i 
and j , then we can estimate the noise correlation function 
R(9) of the 19 GHz data (after removing the monopole, 
dipole, and quadrupole) by R(6) = Ng 1 ^2 ij y i yj/a i (jj, 

where y\ is the antenna temperature of the i th pixel and 
the sum is over all N$ pairs of pixels separated by 9. Sub- 
stituting these relations into eq. (||) gives estimates of Ad 
which are from 6% to 17% larger than those in Table 1. 



3.1.2. Monte Carlo Simulations 

As another test of the robustness of the estimates of d, 
we employed Monte Carlo simulations in which the tem- 
plate maps are sliced into eighteen regions of equal area, 
each corresponding to a range of Galactic latitude \b\. In- 
side each of these regions the pixels are rearranged in ran- 
dom order, so that the latitude dependence is preserved 
but the longitudinal correlations are destroyed. Repeat- 
ing this procedure 1000 times yields distributions of dVs 
consistent with zero mean and with standard deviations 
A&M ~ 1.6 times larger than the formal errors Ad (see 
Table 1). 

a simple noise- weighted least-squares fit. 
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3.1.3. Sky Rotations 

Because of the approximate axial symmetry of the 
Galaxy it is natural to ask if the correlations are sim- 
ply due to overall large-scale Galactic structure common 
to all emission components. To test this hypothesis, we 
repeated the analysis with 2 x 2 x 36 = 144 transformed 
maps, rotated around the Galactic axis by multiples of 10° 
and/or flipped vertically and/or horizontally. 

For a 20° Galactic cut, the correct template map has 
the highest of all 144 correlations, while the distribution 
of these correlations has standard deviation A&t ~ 2.5 
times larger than the formal errors Ad. Likewise, the cor- 
rect DIRBE maps have the highest of all 144 correlations 
and the standard deviations are A&t ~ Acxm, even for 
a 30° Galactic cut (see Table 1). In contrast, we find 
no significant correlation between the 100/im and Haslam 
maps, indicating that synchrotron and dust emission are 
not strongly correlated at high latitudes. 

Since both A&m and Ackt have correlated signal con- 
tributing to the noise, they constitute overestimates of the 
true error bars. 

3.2. Latitude Dependence 

To investigate the dependence of the correlation on 
Galactic latitude, we sliced the maps into six regions of 
equal area, each corresponding to a range of latitude 
Figure shows the results for the 100/xm map. Note that 
d from the 19 GHz DIRBE correlation is almost the same 
at each latitude band, indicating that this correlation is 
not dominated by one or two nearby clouds. The 19 GHz 
Haslam correlation is found to be more concentrated in 
the Galactic plane 2 . 
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Fig. 2. — The dependence of a on Galactic latitude \b\ for the 
100/xm map. The error bars are given by Aaj. 

3.3. Scale Dependence 

One way to determine angular scale of the correla- 
tion is to compute the angular cross-correlation func- 
tion (CCF) of the 19 GHz map with the template 
maps. We computed this according to CCF(9) = 



T l t 1 X Gal, j y l /0-j I Ei,' " 



, where the sums are again 



over all pairs of pixels separated by 9. 

As an example, Figure || shows the CCF(9) of 19 GHz 
and 140^m map, as well as the auto-correlation function 
of 140/im map. It is clear that the correlated structure in 



the two maps is on small angular scales (9 < 10°) and that 
the CCF is well-behaved on all angular scales. The latter 
gives us additional confirmation that there are no large, 
unknown systematics that compromise the analysis. Simi- 
lar results are found for the other two DIRBE maps. The 
lower signal-to-noise for the correlation of the two radio 
templates result in CCF's which are less demonstrative. 
Note that CCF(0) oc d. 
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Fig. 3.— CCF(0) between 19 GHz and 140 map for a 30° 
Galactic cut, together with the 140 fim auto-correlation. 

As another way to investigate the dependence of the 
correlation on the angular scale, we repeated the analysis 
after high-pass filtering the 19 GHz map using a partial sky 
multipole removal technique (see Tegmark & Bunn 1995). 
Only after removing spherical harmonic components with 
I > 8 does the correlation significantly decrease, which 
again indicates that large-scale Galactic structure is not 
the source of the correlation. Both types of filtering there- 
fore indicate that the bulk of the correlations are caused 
by fluctuations on scales of several degrees, which is con- 
sistent with the shape of the CCF in Figure [§[ 

4. CONCLUSIONS 

The two synchrotron templates are found to be 
marginally correlated with the 19 GHz map, while all three 
DIRBE far-infrared templates show a significant correla- 
tion. When a synchrotron and a DIRBE template are si- 
multaneously fit to the 19 GHz map (via linear regression) 
the correlation coefficients do not change significantly from 
those listed in Table 1 and, in addition, the two fit parame- 
ters are essentially uncorrelated (\p\ < 0.05). We conclude 
that the correlations with DIRBE dust emission are in- 
dependent of the correlations with synchrotron emission. 
The amplitude of the signal is much larger than expected 
for ordinary (vibrational) dust emission, as shown in Fig- 
ure U Moreover, there is now good agreement between 
different experiments that this correlated component is 
brighter at lower frequencies. 

So what physical component is this? Two contenders 
have been proposed. K96ab argue on physical grounds 
that free-free emission might be spatially correlated with 
dust. However, the correlations between Ha (which is nor- 
mally a good tracer of free-free emission) and CMB maps, 
and between Ha and the DIRBE maps, are weak (L97; 



2 In order to test how important spatially localized features are for the correlation, we cross-correlated 19 GHz with Haslam using only one 
hemisphere. Most of the correlation was found to come from the northern hemisphere, presumably because of emission features such as the 
North Polar Spur. 
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Frequency (GHz) 

Fig. 4. — Frequency dependence of the fluctuating component 
of the Galactic emission for |fe| > 30°. The DMR correlations are 
represented by filled circles (K96b), Saskatoon data by pentagons 
(dOC97) and 19 GHz by a square with error bar given by A&t 
(this work). Dust emission (vibrational modes only) is normalized 
using the DMR correlation at 90 GHz (K96a). The free-free model 
is normalized to the DMR correlation at 31.5 GHz (the rms free- 
free emission derived from Ha is much smaller than this DIRBE- 
corrclatcd component). The synchrotron emission is normalized to 
the correlation between the 19 GHz map and the 408 MHz Haslam 
map (see Table 1), and is consistent with the 31.5 GHz upper limits 
from K96b. The thickness of each curve corresponds to the normal- 
ization uncertainty. The dashed line is an example of emission from 
rotating dust grains (DL98). All fluctuations have been converted 
to DMR angular scales as ST = eaa\ ,, where a l „ , is the rms of 

° Oat' bai 

each template map and the correction factor is 

e={[Y J (2l + l)if 3 Wp MR ]/[Y J ('2l+l)l f3 Wl}f /2 . 

Here WP is DMR window function, is the window function 
of the experiment to be converted and we assume a fi = —3 power 
spectrum slope. 

Kogut 1997; McCullough 1997). One possibility is the 
presence of an extremely hot (>10 6 K), ionized plasma 
(L97); however, DL98 have argued that this cannot be 
true across the whole sky on energetic grounds. These re- 
sults motivated DL98 to suggest that perhaps it is dust 
after all, but emitting through rotational rather than vi- 
brational excitations. As shown in Figure^, this can give 



a spectrum quite similar to that of free-free emission in 
the relevant frequency range. 

Ha maps with better accuracy are currently being 
made by the Dennison, Gaustad and Reynolds groups and 
should help to settle the issue. Another test is indicated 
by Figure [|: although both the rotating dust and the free- 
free models appear consistent with the available data 3 , 
the former predicts a bump and a downturn around 10- 
20 GHz while the latter predicts a continued rise toward 
lower frequencies. A cross-correlation analysis with lower 
frequency data like the 10 GHz Tenerife map may be able 
to discriminate between these two models. 
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ble Fellowship HF-01084.01-96A from STScI, operated by 
AURA, Inc. under NASA contract NAS5-26555. 

Table 1 — Correlations with the 19 GHz map 
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(a) Correlations for jbj>20° (top) and |6|>30° (bottom). 
(t) a has units fj,K (MJy/sr) -1 for the DIRBE templates, fiK/K 
for the 408 MHz template and /iK/mK for the 1420 MHz tem- 
plate. 

(c) ST= {&±Aa T )a G ai. 
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3 We compute the spectral index between 19 and DMR 53 GHz to be ~ -2.4, and between 19 GHz and 1420 MHz to be ~ -2.8. 



